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ABSTRACT 


In  this  report  we  investigated  the  variation  of  dose  rate  with  solid 
angle  firaction  in  structures  of  simple  geometries  and  compared  the  ex¬ 
perimentally  determined  infinite  field  dose  rate  with  that  calculated 
using  the  formula 

d"  =  B  0.5S  E  [s  (0)  )-S  (oj  ,  )1  +  0-088(1-8  ) 
w  w  I  a'  u'  a'  u  '  w' 

The  geometry  factor  for  skyshine  radiation  S  (co)  was  taken  from  Figure 
B37  of  "Structure  Shielding  Against  Fallout  Radiation  from  Nuclear 
Weapons"  by  L.  V.  Spencer  and  was  applied  to  wall-scattered  radiation 
as  described  in  "An  Engineering  Method  for  Calculating  Protection  Af¬ 
forded  by  Structures  Against  Fallout  Radiation"  by  C.  Eisenhauer.  The 
factors  E,  and  were  taken  from  the  Engineering  Manual 

"Design  and  Review  of  Structures  for  Protection  from  Fallout  Gamma 
Radiation. " 

Except  for  the  thinnest  walls,  the  values  of  experimental  and  calcu¬ 
lated  dose  rates  were  within  20%.  This  was  greater  than  the  estimated 
experimental  error,  and  the  difference  followed  a  pattern  that  suggests 
a  definite  discrepancy  between  experiment  and  calculation.  The  same 
patteim  of  discrepancies  was  found  in  a  direct  comparison  of  experi¬ 
mental  and  theoretical  values  of  S  (u). 
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CHAPTER  1 


INTRODUCTION 

The  prediction  of  dose  rates  to  be  expected  in  a  radiation  shelter  due  to  exter¬ 
nal  sources  presents  a  problem  far  more  complicated  than  that  of  determining  the 
attenuation  of  a  narrow,  collimated  beam  of  radiation  traversing  a  thickness  of 
material  and  being  measured  by  a  collimated  detector.  Tte  dose  rate  in  a  shelter, 
besides  the  narrow  beam  contribution  of  direct  radiation,  will  include  secondaries 
scattered  through  all  angles  and  often  of  comparable  intensity. 

During  recent  years  this  aspect  of  radiation  shielding  has  been  intensively 
investigated  with  particular  emphasis  on  applications  to  shelters.  Of  particular 
importance  has  been  a  series  of  machine  computations  of  gamma  ray  scattering  and 
attenuation  functions  prepared  and  compiled  by  Dr.  L.  V.  Spencer^  using  the  moment 
method  and  Monte  Carlo  techniques  based  on  fundamental  cross-section  data.  Com¬ 
putations  of  the  angular  distribution  of  scattered  radiation  at  various  distances  from 
point  and  plane  sources  yielded  twenty- five  graphs  from  which  flexible,  self- 
consistent  methods  of  analysis  were  derived  by  Mr.  Charles  Eisenhauer  and  Mr. 

L.  N.  FitzSimons  working  in  conjunction  with  Dr.  Spencer.  These  curves,  com¬ 
puted  for  Co-60,  Cs-137,  and  1. 12- hr  fallout  radiation,  have  been  published  in  the 
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so-called  Engineering  Manual  ,  which  is  designed  to  assist  architects  and  engineers 
in  the  design  of  structures  for  fallout  protection. 

To  date ,  these  methods  have  been  experimentally  evaluated  only  by  whole  building 
exposures  —  structures  so  complicated  as  to  yield  little  basic  information,  though 
providing  significant  confirmation  of  the  general  applicability  of  the  method.  The 
basic  premises  of  this  important  field,  therefore,  have  not  been  tested  experimen¬ 
tally.  This  experiment  was  planned  to  test  the  premises  by  using  the  simplest  possi¬ 
ble  structure,  thus  reducing  the  measured  dose  rate  to  a  function  solely  of  the  sub¬ 
tended  solid  angle  and  the  thickness  of  the  scattering  walls. 

A  cylindrical  shell  at  the  center  of  a  circular  field  was  used  to  provide  this  sim¬ 
plicity,  with  the  additional  advantage  of  circular  symmetry,  making  possible  the 
replacement  of  the  plane  source  by  a  successira  of  point  sources  along  one  radius. 
Wall  thickness  was  changed  by  adding  successive  half-cylinders,  and  solid  angle 
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was  changed  by  placing  the  detector  at  different  depths  in  a  hole  beneath  the  cylinder. 
Below-surface  positions  were  necessary  to  eliminate  direct  radiation.  The  infor¬ 
mation  derived  from  this  experiment  about  the  variation  of  dose  rate  with  solid  angle 
for  a  given  cylinder  thickness  could  then  be  compared  with  calculations  based  on  the 
Engineering  Manual  and  Spencer's  Monograph.  In  the  second  part  of  the  experiment, 
the  cylinder  was  replaced  by  a  wall  of  1-in.  thick  steel,  6  ft  by  12  ft,  to  test  the 
possible  influence  of  structure  configuration  <m  the  dose  rate. 

The  experiment  is  described  in  detail  in  Chapter  2.  The  basic  data  are  pre¬ 
sented  in  Chapter  3  and  are  analjrzed  in  Chapter  4.  Chapter  5  contains  the  conclu¬ 
sions  derived  from  the  experiment.  An  experimental  curve  for  S  (o))  is  derived  in 

cl 

Chapter  6  and  is  compared  with  the  theory. 


CHAPTER  2 


DESCRIPTION  OF  EXPERIMENT 

The  primary  purpose  of  this  experiment  ^'as  to  investigate  the  validity  of  the 

assumptions  made  by  Eisenhauer  and  Spencer  in  deriving  the  curve  of  G  (cj),  defined 

s 

as  a  scattering  function  such  that  D**  =  ^ich  is  that 

fraction  of  the  upper  hemisphere  not  subtended  by  the  scattering  surface,  and  is 
that  fraction  of  the  lower  hemisphere  not  subtended  by  the  scattering  surface.  This 
requires  experimental  knowledge  of  the  variation  of  dose  rate  with  solid  angle  with 
minimal  structiiral  or  radiation  complications  and  of  the  e^ect  of  wall  thickness  and 
structure  shape  on  the  dose  rate. 


OUTLINE  OF  EXPERIMENT 

To  achieve  a  direct  dependence  on  solid  angle,  the  detector  was  placed  on  the 
axis  of  a  steel  cylinder  subtending  solid  angle  o)  =  1  -  cos  O ,  where  G  is  the  angle 
betw’een  the  cylinder  axis  and  a  line  from  the  detector  to  the  upper  (u>^)  or  lower 
(o;^, )  edge  of  the  cylinder  (Figtire  1).  The  dose  rate  could  then  be  measured  at 
different  detector  positions  along  the  axis,  expressed  in  terms  of  w,  and  a  com¬ 
parison  made  with  G  (oj).  However,  the  detector  had  to  be  placed  below  the  level 
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of  the  contaminated  plane  to  eliminate  direct  radiation,  giving  an  expression 
D*  =  Of  ^  curve  of  derived  from  Spencer’s 

S^(d,a;)  by  the  relation  G^  =  0. 5  |^1  -  S^(d,w)J.  Since  the  G^fw)  curve  in  the  Engi¬ 
neering  Manual  refers  to  fallout  radiation,  calculations  were  based  on  the  S  (d,ci;) 
curve  in  Spencer’s  monograph  for  Co-60  radiation  using  the  formula 


d”  =  0.5  a 


Possible  dependence  of  G  (cu)  on  wall  thickness  was  investigated  by  varying  the  mass 
thickness  of  the  cylinder;  its  dependence  on  structure  shape  was  determined  by  re¬ 
placing  the  cylinder  with  a  rectangular  wall. 
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Figure  1.  Solid  Angles  Subtended  by  Upper  and  Lower  Limits  of  the 
Cylinder  at  Detector  Position  D 

However,  the  majority  of  the  experiments  were  performed  on  the  cylinder  to 
take  advantage  of  its  symmetry,  which  permitted  a  point  source  to  represent  a  ring 
source  in  the  following  wai\  If  a  source  qS^  is  placed  a  distance  r  from  a  detector, 
the  response  will  be 


If  qS  is  considered  a  linear  radiation  density,  then  the  total  response  at  the  center 
®  2 

of  a  circle  of  such  line  charge  is  2»TqS^/r  ,  since  the  contribution  from  one  ele¬ 
ment  is  equal  to  that  from  any  other.  If  this  circumferential  source  is  thought  to 
represent  an  annular  ring  of  area  2aTdr,  then  the  ring  contributes 


dD  =  qS 

o 
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or 


dP 

dr 


2irrD 

m 


and  the  total  field  can  be  approximated  by  a  summation  of  elements  2TrArD^ ,  being 

represented  by  individual  point  sources  a  distance  Ar  apart  along  one  radial  line. 

Since  2!rrD  =  D  is  the  basic  datum ,  all  measurements  are  expressed  in  this 
m  exp  2  ^ 

form  and  are  normalized  to  a  field  strength  of  1  curie/ft^. 

The  same  philosophy  was  followed  in  simulating  the  field  for  the  rectangular  wall 
except  that  the  response  would  vary  with  the  angle  <f>  between  the  perpendicular  to 
the  plate  and  the  source-detector  line.  The  total  response  would  then  involve  a  sum¬ 
mation  over  radiiis  followed  by  a  siimmation  over  the  angle  4>. 


CHOICE  OF  MATERIAL 
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As  Eisenhauer  points  out,  his  curve  for  G  (to)  is  based  on  skyshine;  i.e. ,  radi- 

s 

ation  scattered  by  an  atmospheric  medium.  Practically,  a  denser  scattering  mate¬ 
rial  was  experimentally  necessary  and  generally  desirable  for  shelter  applications. 
The  significant  choices  were  iron  and  concrete,  both  of  which  are  presumed  to  be 
reasonably  equivalent  to  air  as  scatterers,  on  a  mass  thickness  basis,  since  the 
absorption  coefficients  of  the  lighter  elements  are  quite  close  together  at  the  photon 
energies  of  Co-60  or  1. 12-hr  fallout  (around  1  MeV).  Concrete  would  have  greater 
applicability  in  structural  shielding,  but  iron  is  easier  to  shape  into  the  desired  form 
and,  more  important,  is  much  freer  of  voids  and  inhomogeneities.  Since  Monte 

4 

Carlo  calculations  indicate  there  is  little  actual  difierence  in  the  scattering  response 
of  the  two  materials  (Table  1),  iron  was  chosen  for  the  structural  material. 


EXPERIMENTAL  STRUCTURES 

The  basic  cylinder  was  constructed  of  steel  l/2-in.  thick,  5-ft  high,  and  2-ft  in 
diameter.  The  height  seemed  to  approximate  an  infinite  cylinder  closely  enough, 
since  even  at  the  highest  detector  position  was  quite  close  to  0  (specifically  0. 0188). 
This  assiunption  was  justified  by  varying  the  cylinder  height  and  by  observing  the 
change  of  dose  rate.  The  dose  rate  was  close  to  its  asymptotic  value  at  a  5-ft  height 
(Figure  2).  Thickness  could  be  increased  by  adding  1/2-in.  -thick  half  cylinders. 
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TABLE  1 


ROOF  REDUCTION  FACTORS  COMPUTED  BY  THE  MONTE  CARLO 
METHOD  FOR  A  1.  25  MeV  PLANE  ISOTROPIC  SOURCE  ON  A  BARRIER 
OF  30  g/cin2  =  61.5  psf  EFFECTIVE  MASS  THICKNESS 


03 

Concrete 

Iron 

0. 094 

1 

0.  0073 

0. 0076 

0.  234 

0.  0170 

0. 0176 

0. 357 

0.  0245 

0.  0257 

0. 500 

0.  0316 

0.  0326 

0. 658 

0. 0360 

0.  0383 

0.  826 

0. 0391 

0.  0412 

0.913 

0.  0418 

0.  0424 

LOOO 

0.0419 

0.  0429 

Figure  2.  Variation  of  Detector 
Response  with  Cylinder  Height 
for  Various  Source  Positions 
(h  =  3-in.  depth) 
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shown  in  Figure  3,  to  the  side  of  the  cylin¬ 
der  facing  the  source.  Experiments  were 
conducted  for  four  cases  of  effective  mass 
thickness:  (1)  1/2  in.  or  18.6  psf;  (2) 

1  in.  or  37.  2  psf;  (3)  1-1/2  in.  or 
55. 9  psf;  (4)  2  in.  or  74. 5  psf. 

The  rectangular  wall  was  a  1-in.  - 
thick  steel  slab  12  ft  long  by  6  ft  high.  At 
a  detector  position  on  a  line  4  in.  behind 
the  wall,  these  dimensions  gave  solid 
angle  subtended  by  the  upper  edge  of  the 
same  order  as  that  subtended  in  the  cylin¬ 
der.  The  assiunption  that  this  would  give 
a  response  close  to  that  of  an  infinite  plate 
was  tested  and  verified  by  adding  a  small 
extra  plate  (Figures  4  and  5). 

The  cylinder  v.’as  placed  over  a  hole 
of  the  same  diameter  and  deep  enough  for 

■  «  s  s  «  c 
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Figure  3,  1/2-in.  Thick  Cylinders  and  Semi- Cylinders  Used  in  the  Experiments 


a  4-fl  detector  depth.  The  upper  rim  of  the  hole  was  lined  with  lead  bricks  to 
reduce  edge  effects.  The  rectangular  wall  was  alongside  a  hole  that  only  extended 
3  ft  on  either  side  of  the  detector  position,  thus  blocking  some  of  the  wall  from  the 
detector  line-of-sight.  The  change  this  made  in  the  solid  angle  of  the  wall  at  the 
detector  v\'as  significant  only  tor  the  deepest  three  positions  (Table  2).  The  rim 
was  lined  with  lead  bricks,  and  a  steel  plate  covered  the  portion  of  the  hole  behind 
the  detector  to  reduce  airscatter. 
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Figure  4.  1-in.  Thick  Rectangular  Wall  with  Dosimeters 
in  Position 


Figure  5.  1-in.  Thick  Rectan¬ 

gular  Wall  Showing  Detector 
Pit  and  Ratemeter 


TABLE  2 

CHANGE  IN  SOLID  ANGLE  DUE  TO  LIMITED  DETECTOR  HOLE 


Detector  Depth 

L 

(in.) 

“l 

"2 

"2 

36 

0.  0157 

0. 0181 

0.  867 

24 

0.  0312 

0.  0335 

0.931 

12 

0. 0787 

0.  0805 

0.978 

6 

0. 1635 

0. 1637 

0.999 

12  ft 


=  solid  angle  fraction  subtended  by  the  portion  of  the  plate 
not  hidden  by  the  edges  of  the  detector  pit 

032  =  solid  angde  fraction  subtended  by  the  entire  plate. 
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POINT  SOURCES 


Point  sources  were  placed  manually  at  selected  radial  positions.  Sources  of 
approximately  20,  7,  and  0.5  curies  were  used,  depending  on  detector  depth  and 
source  radius.  The  0. 5  curie  source  was  rod>mounted  and  positioned  by  hand;  the 
7-curie  and  20-curie  sources  were  cranked  into  position  through  polyethylene  tubing 
by  a  cable-type  retrieval  unit  used  for  source  positioning. 

INSTRUMENTATION 

At  the  deeper  positions  with  smaller  dose  rates ,  the  detector  used  was  an  ion- 
izaticm  rate  chamber  removed  from  a  standard  Victoreen  Model  592  Survey  Meter. 

3 

This  chamber,  with  an  approximate  volume  of  320  cm  ,  was  modified  by  increasing 
the  resistance  across  the  chamber  by  a  factor  of  6  to  provide  greater  sensitivity. 
The  chamber  was  enclosed  in  a  polyethylene  cover,  1/8- in.  thick,  to  ensure  that 
electrons  caused  by  free-air  ionization  did  not  penetrate  the  chamber  wall.  The 
output  of  the  chamber  was  fed  to  a  solid-state  amplifier  ~  a  specially  designed 
low-noise,  feedback  type.  To  remove  statistical  fluctuations,  the  output  of  this 
amplifier  was  fed  through  a  time-integrating  circuit  to  a  John  Fluke  digital  null  volt¬ 
meter.  In  this  manner,  the  output  of  the  chamber  could  be  read  to  approximately 
two  parts  in  a  thousand  with  excellent  reproducibility. 

Near  the  surface,  Victoreen  Model  362  pocket  dosimeters  (200  mr)  were  used 
because  of  the  appreciable  vertical  extent  of  the  ratemeter.  These  were  read  with 
a  portable  charger- reader  designed  and  cxxistructed  by  Tech/Ops. 

CALIBRATION 

The  ionizati(xi  rate  chamber  used  in  the  experiment  was  calibrated  by  exposing 
it  to  a  standard  Co-60  source  previously  calibrated  by  the  National  Bureau  of  Stand¬ 
ards.  Both  chamber  and  source  were  placed  on  an  essentially  massless  calibration 
bench  approximately  9-1/4  ft  above,  and  parallel  to,  the  ground.  Voltage  readings 

of  the  chamber-amplifier  output  were  taken  at  various  source-to-detector  distances 

2 

ranging  from  4  ft  to  11  ft.  These  data  were  corrected  for  air  attenuation,  R  atten¬ 
uation,  and  reflection  of  radiation  from  the  ground  to  the  detector.  The  resulting 
value  of  sensitivity  obtained  by  extrapolating  these  calibration  data  to  zero  radius 
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for  the  detector-amplifier  combination  was  273  V  output/r/hr.  The  detector  re¬ 
sponse  was  determined  to  be  essentially  linear  within  1/2  of  1%  for  voltage  below 
10  V.  Thus,  when  taking  data  for  this  series  of  experiments,  we  vary  the  source 
strength  to  keep  the  maximum  detector  output  below  this  value  as  much  as  possible. 

A  sectmdary  calibration  of  the  detector  was  performed  by  placing  it  within  the 
hole  and  locating  a  small  source  of  Co-60  in  a  fixed  position  relative  to  the  detector 
immediately  after  tiie  fundamental  calibration  was  performed.  This  method  of  sec¬ 
ondary  calibration  was  repeated  to  evaluate  the  detector  performance  at  the  start 
and  conclusion  of  each  day's  experimentation. 

The  calibraticm  of  dosimeters  necessary  to  convert  microampere  readings  from 
the  charger-reader  to  dose  values  in  milliroentgens  was  conducted  outdoors  to  re¬ 
duce  scattering  into  the  dosimeters-  The  dosimeters  were  carefully  placed  at  heights 
above  the  ground  that  would  minimize  local  efiects  of  ground  scattering  without  un¬ 
duly  cmnplicating  the  procedure.  All  exposures  were  to  Co-60  sources  of  knoun 
strength.  Groups  of  identical  dosimeters  were  placed  a  given  distance  from  the 
source.  Different  exposures  were  obtained  by  varying  the  time  of  exposure.  The 
dose  was  calculated  based  <m  an  inverse  square  law  behavior  for  source-dosimeter 
distance  and  no  scattering,  since  the  distances  involved  were  small. 

OPERATIONAL  PROCEDURES 

Personnel  engaged  in  the  experiment  were  equipped  with  200  mr  direct- reading 
dosimeters  and  film  baches ,  and  both  the  air-supported  experimentation  building 
and  the  ccmtrol  office  were  continuously  monitored  by  Tech/Ops  Model  492 
"Gammalarms.  *  In  addition,  during  all  manipulation  of  sources,  two  portable  sur¬ 
vey  meters  were  used  for  personnel  monitoring. 

Dosimeters  were  accurately  positioned  at  depths  of  1,  2,  3,  4. 5  ,  and  6  in.  cm 
the  axis  of  the  cylinder.  They  were  clipped  to  a  plastic  holder  and  exposed  simul¬ 
taneously.  While  one  operator  mcmitored  from  a  safe  distance,  the  other  attached 
the  point  source  to  a  14-ft  handling  rod  and  positioned  it,  in  the  case  of  the  0. 5-cmrie 
source ,  or  fastened  the  end  of  the  polyethylene  tube  at  the  prc^r  position  and 
cranked  out  the  larger  sourcjes.  After  the  exposure  was  completed,  the  source  was 
replaced  in  its  storage  container,  and  the  dosimeters  read  and  replace  for  the  next 
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run.  Readings  were  taken  at  radial  positions  of  1. 75,  2. 25,  2. 75,  3.  25,  3.  75, 

4.25  ,  4.75  ,  5.25  ,  5.75  ,  6.5  ,  7.5  ,  8.5,  11.0,  13.0,  15.0,  17.0,  19.0  ,  22.5  ,  27.5, 
32.5,  37.5,  and  42.5  ft,  and  extrapolations  were  made  at  1.0,  9.5,  and  45  ft. 

The  ratemeter  was  mounted  on  an  aluminum  rod  and  positioned  by  a  pulley 
system.  The  same  set  of  radial  positions  was  measiired  at  depths  of  2.  75  ,  4. 25, 
5.75,  8.75,  11.75,  17.75,  23.75,  35.75,  and  47.  75  in.  The  obsciure  depth  settings 
resulted  from  an  initial  error  of  0. 25  in.  in  the  ratemeter  position.  The  solid  angle 
subtended  by  the  detector  varied  with  its  depth  as  u;  =  1  -  cos  0 ,  where  ©  is  given 
by  tan  0  =  a/h;  a  is  the  radius  of  the  cylinder;  and  h  the  detector  depth. 

For  the  rectangular  walls,  dosimeters  were  placed  at  depths  of  1. 0,  2. 0,  4.  0, 
and  6  in.  and  the  ratemeters  at  depths  of  4. 0,  6. 0,  12. 0,  24. 0,  and  36  in.  The 
source  was  placed  at  distances  of  1. 75,  3. 75,  7. 5,  22. 5,  and  42. 5  ft,  each  radial 
line  being  repeated  at  0,  20,  40,  55,  70,  and  80  deg.  The  angular  positions  were 
chosen  <m  the  basis  of  a  crude  estimation  of  flux  change  with  angle.  Solid  angle 
measurements  were  based  on  the  inner  c<^es  of  the  structures. 
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CHAPTER  3 


EXPERIMENTAL  DATA 


Tables  3,  4,  5,  and  6  show  measured  values  of  =  2rrD^  for  effective 
mass  thicknesses  of  18.6,  37.2,  55.9,  and  74.5  psf,  respectively.  Values  at 
P  =  1. 00,  9. 5,  and  45  ft  are  extrapolated  from  curves  of  vs  p.  Table  7  gives 
Dexp  for  the  rectangular  wall.  Dosimeter  readings  are  at  depths  of  1  in. ,  2  in. , 

3  in. ,  and  the  first  listed  at  (4a)  4  in.  ;  ratemcter  readings  are  at  4(b),  6,  12,  24, 
and  36  in.  Table  8  lists  the  value  of  solid  angle  fraction  associated  with  each  depth 

2 

in  the  cylinder  and  wall.  All  values  are  normalized  to  a  field  strength  of  1  curie/ft  . 


TABLE  3 

DOSE  RATE  DATA  DC  S-PT  CTLDJDER 
CX«  -  IS.  6  paf)  Df^xp  OBr/hrVVarie/ft) 


P(ft) 

B 

B 

2.75 

4.25 

B 

6.0 

9.75 

1L75 

17.75 

23.75 

35.75 

47.75 

1.00 

fTT 

6.9 

5.2 

5.7 

4.0 

4.3 

2.79 

3.2 

L87 

L220 

0.695 

0.342 

0.149 

0.0790 

1.  T5 

5.  2 

4.  4 

3.59 

3.65 

3. 59 

2.94 

2.24 

2.29 

L39 

0.923 

0.464 

0.265 

0.111 

0.0667 

2.23 

3.9 

3.4 

2.92 

2.90 

2.92 

2.79 

L70 

1  L71 

L092 

0.741 

0.374 

0.217 

0.0940 

0. 0499 

2.73 

3  2 

2.7 

2.33 

2.3 

2.33 

L90 

L44 

L40 

0. 922 

0.621 

0. 320 

0. 191 

0.0833 

0.0439 

3.23 

2.9 

2.3 

2.02 

L33 

2.02 

1.54 

1.36 

1.19 

0.907 

0.553 

0.290 

0  170 

0.0745 

0.0395 

3.75 

2.4 

2,  0 

L73 

1.65 

1.73 

1.39 

1.09 

1.06 

0.709 

0.490 

0.254 

0. 149 

0.0690 

0.0966 

L2S 

2.1 

1.9 

L55 

1.46 

1.55 

1.21 

0.972 

0.934 

0.643 

0.446 

0.231 

0. 137 

0. 0625 

0.0337 

A.  75 

1.  • 

1.6 

L32 

L31 

L32 

L07 

0.944 

0.951 

0.563 

0.392 

0.204 

0.  125 

0.0674 

0.0317 

5.2S 

L4 

L24 

L22 

L24 

0.996 

0. 799 

0.769 

0.514 

0.360 

0. 199 

0.115 

0.0537 

0.0293 

5,  75 

L5 

L3 

LIO 

1.07 

LIO 

0.973 

0.714 

0.699 

0. 474 

0.326 

0. 174 

0.107 

0.0606 

0.0274 

6.5 

L3 

LI 

0.997 

0.949 

0.997 

0. 732 

0.6U 

0.611 

0.416 

0.279 

0.  US 

0971 

0.0456 

0.0249 

7, 5 

L2 

LO 

0.945 

0.954 

C.  9^ 

0.697 

0.554 

0.549 

0.371 

0.257 

0.139 

0.0967 

0.0409 

0.0222 

S.5 

LI 

0.9 

0.792 

0.756 

0.  792 

0.610 

0.497 

0.494 

0. 332 

0. 125 

0.0773 

0.0369 

0.0203 

9.5 

0.96 

0.9 

0.70 

0.66 

0.70 

0.54 

0.439 

0.42 

0.295 

0.211 

0.117 

0.0702 

0.0339 

0.0190 

ILO 

0.  97 

0.  74 

0.579 

0.619 

0.579 

0. 494 

0. 397 

0.412 

0.<«4 

0. 196 

0. 101 

0.0632 

0.030 

0.0166 

13.0 

0.6S 

0.59 

0.515 

0.493 

0.515 

0.394 

0.327 

0.314 

0.217 

0.  U7 

0.0833 

0.0628 

0.025 

0.0139 

15.0 

0.60 

0.51 

0.435 

0.431 

0.435 

0.349 

0.290 

0.279 

0.193 

0.133 

0.0744 

0.0469 

0.0124 

17.0 

0.53 

0.46 

0.399 

0.395 

0.  399 

0.317 

0.239 

0.253 

0.170 

0.121 

0.0661 

0.0410 

0.0197 

o.oloa 

19.0 

0.47 

0.40 

0.349 

0.336 

0.349 

0.271 

0.225 

0.212 

0.152 

0.108 

0.0588 

0.0365 

0. 0177 

0.00988 

22.5 

0.39 

0.33 

0.301 

0. 272 

0.3013 

0.223 

0. 192 

0.179 

0. 129 

0.0919 

0. 0601 

0.03U 

0.0151 

0.00849 

27.5 

0.31 

0.27 

0.239 

0.226' 

0.239 

0. 179 

0.149 

0.143 

0. 103 

0.0729 

0.0407 

0.0255 

0.0121 

0.00704 

32.5 

0.2S 

0.22 

0.199 

0.  192 

0. 199 

0. 151 

0.129  1 

0.119 

0.0849 

0.0589 

0.0334 

0.0211 

0.0104 

0.00671 

37.5 

0.23 

0.194 

0.169 

0.156 

0.  169 

0.133 

0.1099 

0.1039 

0.0733 

0.0618 

0.0299 

0.0179 

0.00890 

0.00604 

42.5 

0.19 

0. 16 

0.142 

0.135 

0. 142 

0.112 

0. 09141 

0.0899 

0.0649 

0.0467 

0.0263 

0.0163 

0.00784 

0.00443 

45.0 

0.19 

0.15 

0.133 

0.13 

0.133 

0. 106 

^^0^ 

0.084 

0.060 

0.0448 

0.0299 

0.0152 

0.00739 

0.00412 
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DOSE  RATE  DATA  IS  5- FT  CYLINDER 
pCg  =  55. 9  psO  Dejm  {inrAr)/(c-jrie/ft) 
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TABLE? 


RATE  DATA  IN'  6  FT  1 12  FT  RECTANGULAR  WALL 
a  =3T.2pi4  D _ (mrAii)/(carie/R) 


l.Ct 

0.28 

0.246 

0.228 

0.196 

0.0720 

1.7S 

0.211 

0.204 

0.164 

•.ISO 

0.0546 

3.75 

0.118 

0.119 

0.0813 

0.0677 

0.0241 

7.5 

0.0634 

0.0002 

0.04M 

0.0319 

•.•102 

22.5 

O.OIK  1 

0.0144  1 

0.00076 

0.00060 

0.00337 

42.5 

0.0090  1 

0.00887  1 

0.0M77 

0.0*436 

0.00199 

45. 

O.OOK  1 

0.0065  j 

•.00460 

•  .00420 

o.oom 

h  °  2  ia. 

l.«0 

1 

0.175  1 

0.169  1 

•.156 

0.151 

•  .•650 

1.75 

0.134  1 

0.129 

0.111 

0.112 

0.0451 

3.75 

0.0768  1 

0.0817  1 

•.0578 

0.0514 

0.0197 

7.5 

0.0410  j 

0.0412 

0.0289 

0.0235 

O.OOOO* 

22.5 

0.0122  j 

0.0110 

•.00M9 

0.00633 

•.•0262 

42.5 

0.0070 

0.00629  1 

•.00347 

•.0«336 

•.••173 

45. 

0.0068  1 

0.00610 

0.00339 

•.••332 

0.00166 

•.ears 

e.Mos 

e.am 

•.Mmi 

•  .MU7 

•.Ml«7 

0.M104 


•.ou* 

•.•IM 

•.•MW 

•.••m 


1 

1.00  1 

0.1150 

0.105 

1 

0.101  1 

0.104 

o.osw 

0.0120 

1.75 

O.C?SS 

0.0861 

O.OSIS 

0.0634 

0.0378 

O.OUO 

3.75 

0.0528 

0.0537 

0.0403 

0.0387 

0.0169 

O.OIU 

7.5  ! 

0.0281  j 

0.0290 

0.0209 

0.0176 

•  .••656 

O.OOIOO 

22.5 

0.0069 

0.00759  1 

0.00409 

0.00631 

•  .••202 

0.000805 

42.5 

0.0044 

0.00458  1 

•.0026* 

•.••198 

O.OOUO 

•  .•••734 

45. 

0.00425 

0.00450 

0.00252 

0.00182 

•.••116 

•.•••722 

h  3  4a  is. 


1.00 

0.0850 

0.0760 

•.0750 

O.OSOO 

•.•360 

O.OIJO 

1.75 

0.0643 

0.0603 

0.0567 

0.0693 

0.0291 

•  .•120 

3.75 

0.0360 

0.0401 

0.0301 

0.0296 

•.0137 

•.•109 

7.5 

0.0201 

0.0212 

0.0151 

0.0138 

•.••630 

•.••188 

22.5 

•.••72 

0.00591 

0.00396 

•.004« 

O.OOIOO 

0.010696 

42.5 

0.0035 

0.00343 

0.00217 

•.•021S 

O.OOIOO 

0.000560 

45. 

0.0034 

0.00332 

•.00206 

0.00206 

•  .••104 

0.000558 

TABLE  7  (C-oM'd.) 

DtSE  RATE  a\TA  IN  6  FT  x  12  FT  RECTANGl  LAR  WALL 
=  37 . 2  p»f)  <iEr/hr)/(cJiric/fi) 


0  - 

“T— ” 

1 

r - 

0^ 

!  2®“ 

1 

I 

! 

i  55° 

i 

!  70° 

J _ 

90° 

h  =  4b  in. 

1 

1.00  I 

c.o»o 

j  O.O^D 

'  0.0*90 

:  0.0900 

]  0.0595 

0.0125 

1.T5 

0.0734 

0.0752 

:  0.0651 

0.0633 

i  0.0351 

0.0125 

3.75 

0.03§1 

0.02‘i<> 

0.0296 

0.0301 

0.0139 

0.0114 

7 . 5 

0.0202 

i  0.02M 

0.0151 

O.OISS 

0.0054* 

0.00195 

22. D 

O.CK>6SO 

,  0.00577 

j  0.00444 

0.00400 

1  O.OOlCo 

0.00075 

42. S 

0.003&5 

'  0.00363 

■  0.00230 

0.00220 

1  0.00096 

0.00057 

45. 

0.00340 

j  0.00330 

!  0-00221 

0.00212 

i  0-00092 

.i.  _ i 

0.00056 

h  =  6  in. 


1.00 

0.0450 

J-  — 1 

'  0.0^00 

0.0390 

;  0.0460 

i 

0.0250 

0.00900 

1.75 

0.0360 

1  0.0374 

0.0301 

0.0337 

i  0.0201 

0.00900 

3.75 

0.02U 

i  0.0215 

0.0177 

0.0196 

'  0.0101 

0.00999 

7.5 

o.cir: 

:  0.0120 

0.00914 

0.00^53 

0.00356 

0.00122 

22.5 

0.00396 

]  0.003*2 

0.00296 

:  0.00263 

0. 00110 

0.000.504 

42.5 

0.00236 

■  0.00224 

0.00150 

C. 00150 

0.00064 

O.0CO396 

4S. 

0.0021* 

.  0-002U 

0.00145 

0.00145 

1  0.00061 

0.000399 

h  =  12  ir.. 


1.00 

0.0140 

:  0.0130 

,  0.0117 

0.0145  0.0120 

j  0.0105 

1.75 

0.01123 

1  0.0113 

:  0.009*5 

0.0116  ’  0.00936 

!  0.0100 

3.75 

0.0076 

0.00744 

0.0064* 

O.OC72)  0.00524 

,  0.00601 

7.5 

0.00476 

;  0.00497 

0.00376 

0.00370  0.00167 

^  0.000570 

22.5 

0.00179  1 

0.00192 

j  0.000999 

;  0.000922  ,  0.900462 

!  0.0001S2 

42.5 

'  0.00105 

0.00106 

0.000759 

’  0.0006*2  0.000292 

i  0.000163 

45. 

0.00103  i 

0.00103 

0.000750 

1  0.000679  ■  0.009295 

J  i 

j  0.000160 

h  =  24  m. 


l.OO 

I 

\  0.00430 

1  0.00440 

i  0.00340 

;  0.00355 

0.00359 

1  0.0033S 

1.75 

0.00337 

:  0.003W 

■  0.00514 

j  0.00340 

i  0.0037 

1  0.00320 

3.75 

i  0,00307 

I  0.00316 

0 . 002f.5 

;  0.00307 

0.0021% 

^  0.0021% 

7.5 

1  0.002D1 

1  0.00204 

0.00156 

0.O0157 

0.00069 

:  0.00025 

22.5 

0.00032 

1  0.00044 

0.00061 

0-00053 

,  0.00021 

:  0.00010 

42.5 

0.00050  ! 

I  0.00051 

;  0.00035 

0.000315  1 

0.00013 

j  0.000C%5 

45. 

0.00049  : 

0.000499 

!  0.00034 

O.OOG309  i 

_A  ..  J 

0,000125 

1  O.OOOC79 

J  ..  . 

b  =  36  la. 


1.00 

0.00249  j 

'  0.00240 

r -  -f- 

j  0.00138  !  0.00155  | 

0.00125 

1  0.00195 

1.75 

0.00211 

0.00214 

1  0.00135  j  0.00149  j 

0.00126 

!  0.00200 

3.75 

V. 00145  i 

I  0.00150 

i  0.00126  ’  0.00142 

0.00139 

0.0C2D9 

7.5 

0.00110 

0.00110 

1  O.OOC930  I  0.000912  ! 

0.000335  i 

0.000123 

22.5 

0.0004C5  ' 

O.0«!4«5 

{  0.000336  j  0.000292  1 

0-000104 

0.0000517 

42.5 

o.&oo?:2  ! 

0.000283 

1  0.000193  :  0.000157  ! 

0.00;>^552  ^ 

0.0000411 

0.000265  ; 

O.OOOKO 

:  0.000196  j  0.000156  i 

..L.: _  i  ‘ 

0.0000540  ! 

J 

0.0000402 
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TABLE  8 


SOLID  ANGLE  FRACTIONS,  w,  AND  DETECTOR  DEPTHS,  h 
5  FT  CYLINDER  AND  6  FT  x  12  FT  STRUCTURE 


Cylinder 

Rectangular  Structure 

h(m. ) 

a: 

Ul 

/ 

u 

u 

u 

u 

1.00 

0.0188 

m 

0.01426 

0.808 

2.00 

0.01821 

mm 

0.01388 

0.670 

2.75 

0.01777 

0.776 

3.00 

0.01767 

0.757 

0.01356 

0.556 

4.00 

i 

0.01330 

0.465 

4.25 

0.01704 

0.666 

4.50 

0.01686 

0.649 

5.75 

0.01633 

0.568 

6.00 

0.01612 

0.553 

0.01250 

■SSI 

8.75 

0.01497 

0.410 

11.75 

0.01367 

0.286 

12.00 

0.01120 

0.168 

17.75 

0.01165 

0.171 

23.75 

0.01007 

0.107 

24.00 

0.00862 

0.0710 

35.75 

0.00767 

0.052 

36.00 

0.00714 

0.0385 

47.75 

0.00610 

1 _ 

0.030 

CHAPTER  4 


ANALYSIS  OF  EXPERIMENTAL  DATA 
INTRODUCTION 

The  calculational  technique  based  on  Spencer's  data  has  proved  very  useful  in 
evaluating  and  predicting  protection  factors  of  various  shelters  within  the  rather 
wide  experimental  limits  and  generally  under  conditions  requiring  uncertain  exten¬ 
sion  of  the  theor>'.  It  is  therefore  desirable  to  have  an  unambiguous  test  of  the  basic 
data  and  assumptions  that  have  g(Mie  into  the  curves  of  the  Engineering  Manual. 

2 

It  will  first  be  necessary  to  outline  the  methods  of  the  Manual .  and  their  deri¬ 
vation  from  Spencer's  work,  and  then  to  demonstrate  how  the  data  that  have  been 
collected  can  be  compared  with  this  information.  The  method  of  calculating  the 
infinite  field  dose  rate  from  the  experimental  data  and  from  the  theoretical  curves 
is  described.  The  experimental  formulas  are  given  and  applied  to  the  data. 

DETERMINATION  OF  DOSE  RATE  AS  A  FUNCTION  OF  SOLID  ANGLE 

The  radiation  reaching  a  detector  located  in  a  structure  can  be  divided  into 
three  components:  direct,  air-scattered,  and  wall -scattered.  The  direct  radiation 
is  the  straight-line  component  from  source  to  detector;  the  component  of  air-scattered 
radiation  that  has  been  scattered  through  large  angles  is  only  slightly  penetrating 
and  ot  little  consequence  in  a  structure  with  thick  walls  and  limited  open  ceiling;  the 
small-angle  component,  introduced  into  calculations  by  the  buildup  factor,  is  near 
the  primary  in  energy  and  penetration. 

The  wall-scattered  radiation  has  been  deflected,  by  one  or  more  Compton  scat¬ 
terings  ,  to  the  detector  position  and  can  usefully  be  considered  as  originating  at  the 
walls.  Such  an  approach  is  valuable  because  the  dose  recei^^ed  at  a  detector  can  then 
be  considered  proportional  to  the  solid  angle  of  the  walls  viewed  by  the  detector  if 
the  walls  have  a  uniform  distribution  of  radiation.  Such  an  assumption  is  justified 
by  assuming  that  the  differences  in  distribution  of  impinging  radiation,  above  and 
below  the  detector  line,  will  be  masked  by  multiple  scattering  in  the  wall;  therefore, 
equal  solid  angles  subtended  above  and  below  the  detector  line  will  yield  equal  detec¬ 
tor  responses.  The  Engineering  Manual  curves  are  based  on  this  assumption,  which 
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would  be  expected  to  introduce  the  greatest  discrepancy  between  experiment  and 
theory. 


If  this  assumption  is  accepted,  the  dose  rate  at  the  detector  in  a  given  structure 

will  be  a  fimction  only  of  the  solid  angle  subtended  by  the  upper  and  lower  limits  of 

the  structure.  Then  a  scattering  function  G  ((i;)  can  be  postulated  such  that  the  dose 

8 


rate  d"  =  afX^)  )j, 

and  shape  of  the  structure. 


where  a  is  a  function  of  the  mass  thickness 


3 

Eisenhauer  has  argued  that  wall-scattered  radiation  is  close  enough  to  air- 

scattered  radiation  in  angular  distribution  to  be  approximated  reasonably  well  Iqr  the 

already  calculaU^  distribution  function  for  air-scattered  radiation.  Then  the  wall- 

scattering  function  G  (cj)  is  defined  as 

8 


u 


t(3\  cos  6)  d(co6  O) 


-(1-w) 


(1) 


£(X,  cos  0)  is  plotted  on  p.  84  of  Spencer's  monograph.  ^  More  directly,  the  curve 
for  G  (a;)  found  in  the  Engineering  Manual  is  derived  from  that  for  S  (  3',ci;)  in  the 

S  8 

Monograph  by  the  relation 


Ggfcj)  =  0.5  1^1- 

The  structural  constant  a  has  been  broken  down  into  B^(X^)S^(Xg)  E.  A  barrier 
shielding  factor  for  the  attenuation  of  radiation  by  the  wall  is  The  curve  in  the 
Engineering  Manual  is  the  same  as  that  for  W(X,d)  in  the  Mon<^raph  multiplied  by  2. 
A  weighting  factor  is  described  on  p.  32  of  Eisenhauer's  explanatory  treatise 
and  represents  the  fraction  of  radiation  scattered  by  the  wall.  Chart  7  of  the  Engi¬ 
neering  Manual  plots  this  weighting  factor. 

The  shape  factor  E  is  described  by  Eisenhauer  on  p.  31  and  is  plotted  in  Chart 
8  of  the  Manual.  This  factor  varies  from  1  for  a  structure  with  infinitely  long  walls 
to  '^2  for  a  square  structure.  The  variation  can  be  interpreted  as  the  additional 


» 
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field  that  would  be  viewed  by  a  detector  in  a  less  eccentric  structure  through  the  two 
walls  (in  a  i^ctangle)  that  would  be  removed  to  infinity  in  the  structure  with  walls  of 
infinite  length. 

Thus  a  complete  system  of  evaluation  has  been  worked  out  based  on  the  calcu¬ 
lated  distribution  of  radiation  intensities  scattered  into  a  particular  solid  angle  by  a 
scattering  interface.  This  is  expressed  in  the  curve  for  G  (u;)  or  S  (  S'.u)  modified 
by  expressions  for  wall  thickness  and  structure  shape. 

COMPARISON  OF  THEORY  WITH  EXPERIMENT 

If  these  theoretical  calculations  have  basis  in  fact,  experimental  knowledge  of 
the  variation  of  dose  rate  with  solid  angle  should  allow  the  comparison  of  those  values 
with  the  infinite  field  dose  rate  computed  from  the  curve  of  S  ( S'.u).  The  simplest 
method  of  comparison  is  to  insert  the  calculated  values  of  and  into  the  Engi¬ 
neering  Manual  curve  and  to  compare  the  result  with  the  measured  value.  This  gives 
an  estimate  of  the  reliability  of  the  composite  expression 

d"  =  0. 5  S  B  E  fs  (w  ,  )  -  S  (w  )  1 
w  w  1  a'  u  '  a'  u'J 

Direct  testing  of  the  theoretical  values  of  S  (u)  requires  additional  assumptions 

Si 

to  derive  an  experimental  curve.  A  discussion  of  this  direct  testing  can  be  found  in 
Chapter  6  of  this  report. 

CALCULATION  OF  INFINITE  DOSE  RATE 

The  ideal  situation  involves  infinite  structures  surrounded  by  infinite  fields  of 
radiation.  The  field  is  in  actuality  limited  by  time  and  the  bounds  of  the  testing  area. 
Experimental  determination  of  a  rapid  convergence  toward  complete  response  with 
structure  height,  expressible  by  the  approximation  of  the  solid  angle  subtended  by 
the  upper  limit  to  zero,  indicated  that  the  practical  limits  of  the  structure  were 
adequate. 

Tte  problem  of  a  finite  cylinder  can  be  accounted  for  in  the  comparison  with 
S  ( 3',a;),  but  the  contribution  of  the  field  beyond  45  ft  amounted  to  a  third  or  more 
of  the  total  dose  rate  and  required  carefiil  estimation  based  on  theory  and  on  data 
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from  the  field  within  45  ft.  As  has  been  discussed  in  Chapter  2  of  this  report,  a 
plane  circular  field  about  a  cylinder  can  be  represented  by  a  succession  of  radial 
points.  The  measured  dose  rate  from  er^ch  point  is  multiplied  by  2sTAr  to  give  the 
dose  rate  from  an  annular  ring.  The  total  dose  rate  from  a  plane  source  is  then 
found  by  summing  these  rings.  The  only  inaccuracy  lies  in  the  finite  size  of  Ar. 

An  additional  inaccuracy  is  introduced  in  the  experiment  with  the  rectangular  wall 
when  the  field  is  simulated  by  a  succession  of  radial  lines  at  different  an^es  to  the 
wall  and  is  separated  by  a  finite  A<». 

As  was  shown,  the  contribution  from  each  annular  ring  is  equivalent  to  the  deriv¬ 
ative  of  the  total  dose  with  distance  dD/dr  at  that  radial  distance.  Then  a  point 
source  qS^  at  a  distance  r  from  a  detector  on  the  axis  of  a  cylinder  will  represent 
an  area  source  to  give  a  gradient 


dr 


2JrqS^e‘^*’B(Mr)  r 


[gs(“u>  '  ^sK'  > 


(3) 


where  gg(w)  is  a  general  scattering  functicm,  assumed  independent  of  r  for  large  r, 
and  B(pr)  is  a  buildup  factor  generally  dependent  on  r. 

If  a  functional  form  can  be  found  for  B(pr),  Eq.  (3)  can  in  theory  be  integrated 

over  a  range  of  r  tlrough  which  r  is  large  enough  to  have  given  Ag  its  independence 

8 

of  r.  Then  the  contribution  from  this  area  to  the  dose  rate  will  be 


2TqS  B(Mr) 

- - -  e^*^dr  •  Ag^ 


(4) 


If  it  is  found  that  Ag^  reaches  its  asymptotic  value  at  a  distance  equal  to,  or  less 
than,  the  greatest  experimental  distance  r  ,  then  r,  can  be  set  equal  to  r  ,r„  =  ««> 

O  X  O  M 

and 


D 


OO 


=  D(r  =  r^)  +  2)rqS^ 


qo 

I 


B(ur)dr  •  Ag^ 


r 


o 


(5) 
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wl^re  D(r  =  r  )  is  the  dose  rate  from  a  field  out  to  r  .  Thus  the  infinite  field  dose 
o  o 

rate  can  be  calculated  if  BQir)  is  known  analytically  and  is  reasonably  simple  and 

if  the  asjnmptotic  value  of  Ag  can  be  found. 

8 

In  traversing  a  medium,  in  this  case  air,  a  beam  of  radiation  is  attenuated  by 
two  principal  mechanisms  at  this  energy:  (1)  photoelectric  absorption,  in  which  the 
phoUm  is  permanently  removed  from  the  beam,  and  (2)  Compton  scattering,  which 
can  occur  a  number  of  times  and  which  may  eventually  scatter  the  photon  back  into 
the  solid  ang^e  viewed  by  the  detector.  In  a  narrow-beam  attenuation  experiment, 
the  total  absorption  coefficient  for  a  particular  material  and  photon  energy  is  meas- 
iu*ed,  but  this  coefficient  in  fact  is  made  up  of  an  energy  absorption  and  a  scattering 
coefficient.  The  energy  absorption  process  follows  a  straight  exponential  decrease 
for  each  element  of  the  radiation  field,  but  this  must  be  modified  by  a  factor  repre¬ 
senting  the  radiation  scattered  back  into  the  beam,  the  buildup  factor.  Thus, 


D  =  D  e 
o 


B(Mr) 

2  ' 
r 


5 

Experiments  on  the  air-scatter  of  Co-60  radiation  have  led  to  a  least-squares, 

0  72 

best-fit  linear  approximation  of  the  buildup  factor  B(pr)  =  1+0. 7(pr)  '  .  This 

was  chosen  as  one  expression  for  B(^r)  with  p  =p^.  As  another  extreme,  B(pr)  =  1 
was  chosen,  assuming  that  no  scattered  radiation  rejoined  the  beam;  a  third  choice 
was  B(pr)  =  1,  P  =  P . ,  assuming  no  scattering  took  place  at  all.  It  was  hoped  that 

OO 

the  experimental  values  of  D  found  by  each  of  these  approaches  would  be  close 
together  and  would  bracket  the  correct  expressicm  among  them.  The  three  approaches 
then  were: 

0  72 

1.  B(pr)  =  1+0. 7(pr)  '  p  =  p^  =  total  (narrow  beam)  absorption  coefficient 

2.  B(pr)  =  1  P  =  P j 

3.  Bdpr)  =  1  p  =  p^  =  energy  absorption  coefficient. 

As  a  derivation  of  asymptotic  scattering  function  g  (cu),  Eq.  (5)  can  thus  be  integrated 

8 

and  used  if  Ag^  can  be  calculated  from  the  data  available  and  is  found  essentially  in¬ 
dependent  of  r  at  the  furthest  measurement  made. 
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The  dose  rate  at  a  detector  <m  the  cylinder  axis  from  wall-scattered  radiation 
may  be  written: 


u 


u 


2t 


dr 


an^o  I  <*“  I 

Cl)  ,  0 


(6) 


where 


o 

q 


=  radial  distance  from  axis  to  source 

=  RHF  =  14. 0  r/hr/curie  at  1  ft  from  Co-60 

2 

=  source  density,  curies/ft 

.2 


2xqS  =  87. 92  for  q  =  1  curie/ft 


o 

H  = 

P  = 


.-1 


absorpticm  coefficient,  ft 
slant  distance  from  the  scHirce  to  the  scattering  area,  ft 


2  2  a^  a  T 

r^  +  h^  + - g - 2  (h  cos  0  +  r  sin  O  cos  ^ 

sin^  0 


> 


e 

h 


a  = 


cylinder  height 
detector  depth 
cylinder  radius 


(i) 


=  1  - 


=  1- 


dose  scattering  function,  steradian 
h  +  e 


-1 


J(b  +  e)^  + 

h 

+  a^ 


•  U 


I 


c 
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As  r  becomes  much  larger  than  a  and  p  r,  independent  of  O  and  then 
Eq.  (6)  becomes 


dD  ZarqS^e-^'-BOir)  >  ^ 


3F 


(7) 


V  ® 


The  integral  in  Eq.  (7)  is  independent  of  r.  If  a  function  g  (X  .w)  is  defined 

8  8 


such  that 


01  2r 


gg(Xe,u)  =  J  J 


(8) 


then 


®  ^  2TqS„e-**'^B«.r)  p 
dr  r 


88<*e-"o>-8s<*e'“y'> 
L  J 


(9) 


where  g^pf^.o;)  represents  the  scattering  properties  of  the  part  of  the  cylinder 
above  cu,  in  a  field  of  radiation  b^inning  a  distance  r^  firom  the  cylinder  axis. 

The  quantity  actually  measured  in  this  experiment  is  dD/dr.  If  the  assumptions 
made  above  are  correct,  the  experimental  data  for  a  g^ven  detector  position  multi- 

{died  by  r/  |^2rqS^e'^^B0^r)J  should  asymptotically  approach  a  constant  value  for 

radii  large  with  respect  to  the  cylinder  radius.  This  value  is  Ag  ,  which  can  then 

8 

be  used  in  Eq.  (5)  to  find  the  infinite  field  dose.  The  result  of  this  computati<»  for 
the  1/2-in.  thick  cylinder  is  shown  in  Figure  6.  The  curves  iire  found  to  approach 
an  asjrmptotic  value  as  r  increases. 

Table  9  lists  Ag  and  D**  for  two  cylinder  wall  thicknesses,  showing  fairly  close 

D 

agreement  among  the  three  approaches  and  the  expected  bracketing  of  the  dose  rate 
incorporating  buildup  by  the  two  extreme  cases. 
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Figure  6. 
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so  35  40  45 


VS  r,  =  18. 6  psf,  BQir)  =  1,  M 


TABLE  9 


ASYMPTOTIC  SCATTERING  FUNCTION  Agg(Xe,aj)  AND  INFINITE- FIELD 
DOSE  RATE  FOR  CYLINDER  WALL  THICKNESSES  OF  18.6  AND  74.5  PSF 
USING  DIFFERENT  FUNCTIONS  FOR  B(ur)  AT  DETECTOR  DEPTH  OF  11.  75  IN. 


Scattering 

Function 

and 

Dose  Rate 

Xg(ps£) 

B(pr)  =  1+0. 7(Mr)®* 

B^ir)  =  1 

M 

BQir)  =  1 
^=^A 

^^s 

ao 

18.6 

0. 0214 

0. 0239 

0.  0230 

D 

11. 17 

10.  26 

11. 54 

^^s 

74.5 

0.  00959 

0.0109 

0.  0105 

d“ 

5.13 

4.  75 

5.33 

•  0«LI»6TOII  •  ■«SS«CHU*CTTt  27 


CALCULATION  OF  D*  BY  THE  ENGINEERING  MANUAL  METHOD 
The  basic  equaticm  accounting  for  wall-scattered  radiation  is 


§-  =  B  S  E}G(a;)-G  (u. ) 
D  w  w  s'  u'  s'  u  ' 


or 


8-  =  0.5B  S  E  S(a),)-S(w) 

D  w  w  a'  u  '  a'  u' 

o  L 

However,  the  experimental  measurements  include  air-scattered  radiation  expressed 
as  the  buildup  factor  in  the  derivation  of  the  far-field  dose  rate  that  must  be  taken 
into  account  as 

D  r  1  0. 088B  r  "1 

#  =  B,(l  -  S^)  )J - 5—  (1  -  S„)  [^S^(a,„, )  - 

(10) 

since  there  is  no  ceiling  scatter  contribution.  The  complete  equati<»  then  becomes 


^  =  0.5B  S  E  [s  (w  ,  )  -  S  (w  )1 
w  w  a'  u  '  a'  u'J 

+  0.088B  (1-  S  ifs  (w  ,)-  S  (w  )|  +  0.088  w  . 

w'  w' t  a'  u  '  a'  u'J  u 


(11) 


The  last  term  represents  the  radiation  scattered  directly  through  the  top  opening  of 
the  cylinder.  Since  it  is  very  small ,  less  than  3%  of  the  total  even  at  the  highest 
position  in  the  thickest  cylinder,  it  can  be  neglected. 

Equation  (11)  implies  quite  a  difference  in  approach  between  the  experimental 
and  theoretical  calculations  for  air-scatter  effects.  Experimentally,  air-scattered 
radiation  was  handled  as  a  multiplicative  factor  modifying  the  radiation  flux  at  the 
scattering  wall ,  ignoring  air-scattered  radiation  striking  the  detector  directly.  The 
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Engineering  Manual  considers  this  radiation  as  an  additive  term  dealing  with  direct 
air-scattered  radiation  <Mily,  attenuated  but  not  scattered  by  the  intervenii^  wall. 
Since  it  would  be  difficult  to  separate  experimentally  the  effects  of  air-scattered  and 
direct  wall- scattered  radiation,  the  Ei^neering  Manual  approach  doesn't  exactly 
mirror  the  physical  situation.  Theoretical  calculations  indicate  that  skyshine  is 
about  10%  of  the  total  for  the  thinnest  (1/2  in. )  walls;  experiments  iialicate  6%. 

There  is  closer  agreement  for  thicker  structures,  pointing  to  a  maximum  error  of 
4%  through  accounting  for  skyshine. 

Equation  (11)  can  be  reduced  to 

§7  =  ®w  >  ■  ^a<“u>]  [o-  5  E  +  0. 088(1  -  S^) j  (12) 

for  final  calculations. 


PRESENTATION  OF  EXPERIMENTAL  FORMULAS 

Dose  rate  was  measured  out  to  42. 5  ft,  and  an  extrapolated  value  at  45  ft  was 

used.  As  shown  in  Figure  6,  Ag  (X  .u)  appears  to  have  reached  a  constant  value 

8  6 

at  this  distance;  also,  the  greater  part  of  D  (~  2/3  of  it)  is  contributed  by  the  radi¬ 
ation  within  45  ft.  This  is  particularly  important  as  the  greatest  source  of  error 

seems  to  be  in  the  determination  of  Ag  . 

s 

Another  major  problem,  that  of  buildup  factor,  was  handled  by  attempting  to 
bracket  it  by  three  sets  of  assumptions: 


1. 

B(pr)  =  1+0. 7(pr) 

2. 

B(pr)  =  1 

3. 

'"ab,  =947.2  11 

B(pr)  =  1 

where  X  is  the  mean  free  path  in  air. 
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This  leads  to  three  separate  solutions  of  the  integral  in  Eq.  (3): 


L  45 


-0.28 

e  ^  dr  =  220. 7 


I  ^ 


-r/445 


Idr  =  2.qS„(^^)  =  159. 


n  -r/947.2  . 

^*•*0  J  - i - 1*  =  ='^‘^0(9^)  = 


and  three  expressions  for  D* 


d"  =  D(P  =  45)  *  220. 7  [^gl(u.„)  -  ) j 

D”  =  D(P  =  45)  -  159. 4  [^g^(u.„)  -  g2(„^  )1 
d"  =  DO.  =  45)  +  221. 3  [^g^(a,^)  -  g3,„^,  ,1. 


finally,  Eq.  (6)  must  be  put  into  three  forms: 


>1  - ; — r-^ - 

■'  87. 92  e"-/"®  [^1  .  0. 7/44S((.r)»- 

[«8<"u>  -  «8<“u'  >]  =  87  32 '-r/445  ^ 

[*»*““’  ■  ’]  '  87. 92  e-*-/®*’- 2  ^ 


to  fiml  the  proper  value  of  g  (w,r  —  «o). 

8 
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Since  the  response  of  the  wall  wotild  be  symmetric  about  the  perpendicular  to  the 
wall ,  measurements  were  made  only  in  one  quadrant  and  the  resulting  dose  rate 
multiplied  by  four  to  simulate  a  structure  with  two  infinitely  long,  parallel  walls. 

The  solid  angle  subtended  by  this  structure  was  computed  using  the  formula  pre¬ 
sented  on  p.  68  of  Spencer’s  Monograph  for  one  surface  and  multipl3ring  this  by  two. 
The  effect  of  the  smaller  size  of  the  detector  bole  or  the  solid  angle  was  found  by 
treating  the  resultant  complicated  shape  as  a  composite  of  simpler  ones. 

EXPERIMENTAL  AND  CALCULATED  VALUES  OF  d“ 

The  raw  data,  presented  in  Table  4.  In  Table  10  D(r)  is  multiplied 

by  Ar  and  is  summed  from  1  ft  to  45  ft  to  gpve  detector  response  firom  a  plane  source 
in  that  area;i.e, ,  D(r  =  45).  In  Table  11  D^^  •  r/2rqS^e”^*^B(lir)  is  presented  for 
the  1/2  in. -thick  cylinder,  and  is  plotted  against  r  in  Figure  6.  Tables  12  and  13 
present  the  asymptotic  values  of  Ag  ..  D(r  =  45),  and  the  values  of  d“*(ci})  corre- 
sptmding  to  the  three  Eq.  (15),  for  each  value  of  X  and  structure  shape.  Calculated 
values  of  D  /D^  are  presented  in  Table  14  and  are  compared  with  the  experimental 
results  in  Table  15. 

The  variation  of  D*’(aj)  with  w  in  the  1  in.  cylinder  is  shown  in  Figure  7.  An 

oo 

extrapolated  value  of  D  at  u  =  1  can  be  derived  from  this  plot  to  give  an  experi- 

r 

mental  estimate  of  B  I  0. 

w  [_ 

method  appeared  to  be  the  approximation  of  the  dose-rate  points  by  an  analytic  func¬ 
tion  of  as  described  on  p.  45.  The  variation  of  dose  rate  with  solid  angle  pre¬ 
dicted  by  the  theory  is  compared  with  the  experimental  values  in  Figure  8. 


5S  E  +  0.088(1  -  S  ) 
w  '  w' 


however,  a  more  consistent 


TABLE  10 

D(r»45).  SLT04ATION  OF  DOSE  RATES  FROM  ANKULAR  KDKJS 
OF  RADIATION  FROM  1  TO  45  FT;  1/2  Di,  CYLINDER 
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TAMLK  U 

ASYMPniTIC  SCATTEIUNG  KrNCTK*N  ASD  INKIMTE  FIELD  Dlisr  KATE 
CYUNDKICAL  STIU  CTl TIES 


lV;;:h 

(in.) 

D<r  =  A5) 

»3' 

u 

1/2-iti.  CvltcMier 

1-00 

33.75 

0.0894 

[  53.48 

[  O.102 

50,01 

0.095 

54.77 

0.917 

2.00 

28.69 

0,07CC 

1  0.088 

42.72 

0.084 

47.50 

0.836 

2.75 

24.44 

0.0682 

1  39. A9 

[  0.078 

36.87 

0.073 

j  40.59 

0.776 

3.00 

23. »S 

0.0638 

38,06 

j  0.072 

35.46 

i  0.068 

j  39.03 

0.757 

A. 25 

19.16 

:  o.o.iAO 

1 

1  0.0613 

28.93 

j  0.059 

5 

1  32.22 

0.666 

4.5 

19.27 

!  0.0524 

1  30. S3 

1  0.060 

28.83 

0.059 

i  32.11 

0.649 

5.75 

13.30 

0.0A3A 

2A.S9 

!  0.050A 

23.33 

i  0.047 

1  25.70 

0.568 

6.00 

15.15 

0.0417 

2A.3-*> 

0.0475 

22.72 

O.OAA 

1  24.89 

0.553 

9.75 

10.13 

0.0299 

16.73 

0.0335 

15.47 

0.0325 

17.32 

0.410 

M.75 

7.03 

0.0214 

11.75 

0.0239 

10.84 

0.0230 

1  12.12 

0.286 

17.75 

3.76 

0.0119 

6.39 

0,0136 

5.92 

0.01.10 

6.6A 

0.171 

23.75 

2.2S 

0.007 AO 

3.91 

0.00939 

3.62 

0.00810 

4.07 

o.io; 

35.75 

1.057 

0.00339 

1.95 

O.OOAOA 

1.699 

0.00395 

1.929 

0.052 

A7.75 

0.579 

0.00202 

1.025 

0.00229 

0.944 

0.00225 

1.074 

0.030 

-in.  Cylinder 

1.00 

33.61 

0.0937 

5A.29 

0.0990 

49.23 

0.09A 

5A.A1 

2.00 

27.67 

0.0761 

AA.A7 

0.0950 

A1.22 

0.091 

A5.60 

2.73 

26.05 

0.069S 

A1.A5 

0.0776 

39. A2 

0.075 

42.65 

3.00 

23. AS 

0.0639 

37.59 
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3A.16 
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A. 25 

20.73 
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33.20 
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34.45 

4.5 
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5.75 
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0.0438 

25.59 
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23.62 
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6.00 
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21.69 

o.OAa> 

23. 6A 

9,75 

9. 963 
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0.0311 

1A.9A 

0.0300 

16.62 

11.75 
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0.0199 

10.70 

0.0212 
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0.0205 

11.06 

17.75 

3.281 

0.0108 
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0.01205  * 

5.202 
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5.892 

23.75 

1.883 

0.00612 

3.23A 

0.00692 

2.986 

0.0065 

3.321 

35.75 

0.939 

0..  287 

1.A72 

0.00316  : 

1.3S3 

0.00305 

1.514 

•1 

0.A53 

0.0015A 

0.7929 

0.00173  j 

0.7288 
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0.9292 

l-l/2-in.  Cylinder 

1.00 

29.67 
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46.20 
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42.74 

0.087 

49.92 

2.00 

2A.S2 

0.0632 

39,47 

0,072  j 

36.00 
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39.35 

2.75 

21.76 

0.0565 

34.23 

0.0610  1 

31.48 
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34.60 

3.00 

21.09 

0.0562 

33.49 
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30.66 

0.057 

33.70 

A.XS 
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25.91 

0.049  \ 

2A.05 
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4,5 
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26.38 
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0.0336 
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20.21  i 

6,00 

13.02 

0.0352 
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0.037  ! 
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20.54 

8.75 

7.623 

0.0210 

12.36 

0.023  1 

11.29 
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12.49  1 

11.75 
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0.0138 
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O.OISS 
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17.75 
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3.91A 
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0.0073 

3.927 

23.75 
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2.299 
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1.013 
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1.065 
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0.5419 
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ASYMI^TOTIC  SCATTKHING  FI  NCTION  ANDISKINH  F  HELD  DDSE  UaTE 
CYUNimiCAL  STIU'CTI  RES 
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0.036 
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19.41 
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0.0293 

i  18.29 

0.033 
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0.032 
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5. 73 
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0.0246 

i  14.54 

0.020 

13.57  ! 

0.029 

15.30 

6.00 

9.454 

e,0253 

■  13.04 

0.029 

i  13.92  j 

0.0-265 

15.32 

5.214 

0.0140 

;  8.502 

0.017 

1  j 

D.016 

9.755 

11.75 

3  2^** 

0 , 0C039 

3.339 

0.0109 

i  4.959  j 

0.0105 

S.S46 

i:.T3 

1.506 

0.00473 

1  2.554 

0.0034 

1  2.367  1 

0.0052 

2.657 

23.75 

0 .  S523 

0.00274 

1.457 

0,0031 

i  1 

0.00295 

1.505 

3s3 , 7  5 

0.3722 

0.00126 

1  0.6503 

0.0014 

j  0.5954  i 

0.00135 

0.6710 

0.2141 
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i  0.387S 
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1. 
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TABLE  13 

INFINITE  FIELD  DOSE  RATES 
RECTANGULAR  STRUCTURE 


Depth 

(in.) 

OO 

1.0 

36.98 

2.0 

26.12 

3.0 

18.52 

4.0 

14.65 

6.0 

9.27 

12.0 

3.803 

24.0 

1.716 

36.0 

0.918 
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TABLE  14 

CALGUL.\TrOKS  USING  THE  ENGINEERING  MANUAL  METHOD 
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TABLE  15 

COMPARISON  OF  EXPERIMENTAL  AND  THEORETICAL  VALLES  OF  D 
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Figure  7.  D“(aj)  vs  oj,  1-in.  Thick  Cylinder 
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Figure  8.  D*  vs  w.  Calculated  and 
Experimental  Values 
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CHAPTER  5 


RESULTS  AND  CONCLUSIONS 
EXPERIMENTAL  ERROR 

The  opportunities  for  error  lay  in  the  interpretation  and  extrapolation  of  the  data 
as  well  as  in  the  experimental  limitations  of  accuracy.  In  taking  readings,  error 
could  arise  from 

1.  Incorrect  source  position 

2.  Surface  irregularities 

3.  Incorrect  detector  position 

4.  Consistent  inaccuracy  of  the  detectors 

5.  Normal  dosimeter  spread  and  reader-charger  readout  limitations 

6.  Finite  size  of  detectors. 

Calculational  errors  would  be  most  likely  in  determining  the  value  of  Ag  from  the 

S 

plots  (see  Figure  6)  and  in  extrapolating  D  (w)  to  w  =  1. 

Most  of  the  experimental  errors  can  be  considered  more  or  less  statistical  fluc¬ 
tuations  because  of  the  large  number  of  independent  positionings  and  readings  tliat 
went  into  the  final  values  of  D*”.  For  the  dosimeters  this  is  true  of  1,  3,  and  5  above; 
for  the  ratemeter,  which  was  left  undisturbed  while  a  radial  series  of  points  was 
measured,  3  could  have  been  serious  if  accurate  positioning  of  the  ratemeter  had  not 
been  ensured.  Surface  irregularities  were  recognized  and  compensated  for.  Since 
the  same  dosimeters  were  consistently  used  in  the  same  positions,  sizable  differences 
in  individual  responses  would  distort  the  final  result.  However,  these  differences 
did  not  exceed  2  to  3%  in  subsequent  measurements.  The  ratemeter,  for  which  6 
might  be  significant,  was  used  only  at  the  greater  depths  where  the  flux  changed  rea¬ 
sonably  slowly. 

Scatter  of  individual  measurements  showed  up  strongly  in  the  plot  of  Ag  vs 

s 

radius  in  Figure  6  to  give  a  somewhat  uncertain  asymptotic  value.  Hopefully,  the 

use  of  three  approaches,  calculated  from  the  Ag  curves  on  different  days  in  different 
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moods,  leads  to  some  statistical  reduction  of  error  in  the  average  value.  The  vari- 

oo 

ation  of  the  quadratic  fitted  to  the  D  (uj)  p>oints  from  those  points  was  no  more  than 
at  the  larger  values  of  oj.  There  could  be  tip  to  4%  error  in  the  treatment  of 
air-scattered  radiation.  On  the  basis  of  these  considerations,  the  data  should  be 
accurate  to  within  10%. 

COMPARISON  OF  EXPERIMENTAL  AND  THEORETICAL  RESULTS 

As  shown  in  Table  15  there  is  good  agreement  between  theoi'y  and  experiment 
for  all  but  the  thinnest  w-all  and  largest  solid  angle.  In  thin  walls,  of  course,  the 
assumption  breaks  down  that  multiple  scattering  will  mask  the  difference  in  response 
to  various  parts  of  the  wall. 

The  G  (cj)  curve  rests  on  the  assumption  that  equal  increments  of  solid  angle 
s 

will  give  equal  response  no  matter  what  the  height  of  the  angles  above  the  source 
plane  or  equivalently  the  detector  position  below  the  plane.  But  in  reality  the 
expression 


oo  '' 

D  =  0. 5  a  1  S  (oj  , )  “  S  (a;  ) 
la'  u  '  a'  u' 

is  valid  only  for  a  detector  position  at  the  level  of  the  source  plane  in  a  semi-infinite 
cylinder;  i.  e. , 


=  1)  =  0.  5  Of  1^1  -  0  =  0.  5  a. 

For  positions  below  the  surface,  the  correct  expression  would  be 

)  =  0. 5  f(d)o  j  .  (16) 

with  f(d)  >  1,  since  lowering  the  detector  is  equivalent  to  raising  the  source  plane 
thus  putting  it  closer  to  each  incremental  scattering  area  that  is  still  exposed.  This 
form,  of  course,  is  necessitated  by  the  retention  of  the  theoretical  S  (w)  rather  than 
by  adopting  a  new,  experimental  form.  The  theoretical  value  would  be  expected  to 
have  sufficient  generality  to  justify  modification  rather  than  replacement.  In  Chapter 
6  the  experimental  curves  are  found  to  tend  toward  Spencer’s  as  a  limiting  value. 
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This  leads  to  two  conclusions  about  the  terms  of  Eq.  (11): 

1.  The  values  of  the  components  of  0.  5S  B  E+0.  088 B  (1  -  S  ) 

I  iV  iV  'n  'n 

may  be  wrong  since  the  theoretical  expression  shovild  be  quite 
accurate  for  ~  1.  Then  if  the  composite  value  is  changed 
to  agree  wth  an  extrapolated  experimental  value  of  D”ataj=  1, 
a  new  set  of  calculated  values  for  D(a;)  is  obtained  as  shown 
in  Table  16.  A  comparison  of  ^"ith  the  theoretical  expres¬ 
sion  is  shown  in  Table  17. 


2.  For  each  measured  value  of  a>,  f(d)  can  then  be  determined, 
assuming  S  (a),d)  actually  represents  wall-scattered  radiation 
and  the  wall-scattered  and  air-scattered  (0. 5  and  0. 088)  com¬ 
ponents  vary  identically  with  depth.  Table  16  also  shows  f(d). 


CONCLUSIONS 


Agreement  within  about  20%  is  fotmd  for  most  of  the  positions  and  wall  thick¬ 
nesses,  except  the  thinnest,  for  which  the  greatest  discrepancy  would  be  expected. 


If  the  validity  of  S  (d.oj)  is  accepted  and  the  barrier  factors  are  questioned,  a 


different  value  for  |  0. 5  S  B  ,  E  +  0.  088  B  (1  -  S  )  |  is  found  from  the  known  be- 


> 


havior  of  S  at  oj  =  1  and  the  extrapolated  value  of  D*.  The  factor  S  (w)  is  defined 
a  “ 

as  equal  to  1  at  w  =  1  and  as  equal  to  0  at  w  =  0.  Therefore,  theory  requires  that  a 
detector  at  the  intersection  of  a  semi-infinite  cylinder  with  a  plane  of  radiation,  sub¬ 
tending  lower  solid  angle  =  1  and  upper  angle  =  0,  receives  dose  rate 


D 


=  or  j^S^(l)  -  S^(0)  =  or. 


Thus  an  experimental  value  for  B  j  0. 5 

L 

for  the  finite  size  of  the  cylinder  and  detector.  The  missing  portion  of  the  cylinder 
between  5  ft  and  infinity  represents  a  very  small  part  of  the  dose  rate  at  the  higher 
detector  positions  and  can  be  adequately  approximated  by  its  theoretical  value. 


S  E 
w 


0.  088(1  -  8^) 


could  be  found,  except 
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TABLE  16 

D®  CALCULATED  FROM  EXTRAPOLATED  VALUES  OF  D*  TO  DETERMINE  f(d)* 

exp  '  ' 


"u' 

D® 

“’calc 

D 

o 

D® 

exp 

D 

o 

f<d) 

D® 

calc 

D 

o 

D 

_j^ 

D 

f(d) 

1/2- in.  C>'linder^ 

2 

1-in.  Cylinder 

0.917 

0. 0947 

0. 1076 

1.  136 

0. 0943 

0. 1092 

1. 158 

0.  836 

0.  0751 

0.  0917 

1.221 

0. 0748 

0.  0895 

1.  196 

0.776 

0.  0650 

0. 0795 

1.  223 

0.0647 

0.  0834 

1.289 

0.  757 

0. 0598 

0. 0766 

1.  281 

0. 0595 

0. 0756 

1.270 

0. 666 

0.0482 

0.0628 

1.303 

0. 0480 

0.0668 

1.392 

0. 649 

0. 0467 

0.0620 

0.0465 

0. 0594 

1.278 

0. 568 

0.  0373 

0.0501 

1.343 

0.  0371 

0.0515 

1.388 

0. 553 

0. 0357 

0. 0490 

1.373 

0. 0356 

0.0483 

1.358 

0.410 

0. 0228 

0. 0337 

1.478 

0. 0227 

0.  0323 

1.424 

0.286 

0. 0149 

0. 0236 

1.584 

0.  0148 

0. 0215 

1.450 

0. 171 

0. 00752 

0.01286 

1.710 

0. 00749 

0.0114 

1.522 

0. 107 

0. 00407 

0. 00787 

1.  934 

0. 00406 

0.  00651 

1.604 

0.052 

0.  00175 

0. 00372 

2.  126 

0.00174 

0.  00296 

1. 697 

0. 000856 

0. 00206 

2.407 

0. 000853 

0.00160 

1.  876 

W  D%).B,[9.SS^E  •  0.9»S<I  -  V*«»  =“  V-u» 

a _ p’m 

1  -  S  (M  )  I  -  0.  OteS  0. 

■  UXD 


1.  D”(1)  -  «2.  M 

o  =  62.40 

=  e. 

2.  D“a).61.»7 

0-62.23 

^  =  0.  IK2 
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TABLE  16  (Cont*d.) 

D*  CALCULATED  FROM  EXTRAPOLATED  VALUES  OF  D  TO  DETERMINE  f(d) 

exp  '  ' 


“u' 

oo 

^calc 

D 

o 

D* 

exp 

D 

o 

f(d) 

O 

o 

D* 

exp 

D 

o 

f(d) 

1-1/2- in.  Cylinder^ 

2-in.  Cylinder 

4 

0.917 

0.  0812 

0. 0930 

1. 146 

0.  0580 

0.0675 

1.  164 

0.  836 

0.  0644 

0. 0774 

1.  203 

0. 0459 

0.  0573 

1.  248 

0.776 

0.  0558 

0.0689 

1.236 

0. 0398 

0.0510 

1.  283 

0.  757 

0.  0513 

0.0674 

1.315 

0. 0366 

0.  0484 

1.321 

0. 666 

0.0413 

0.0521 

1.  262 

0.  0295 

0. 0380 

1. 286 

0.649 

0.  0401 

0.0525 

1.311 

0.  0286 

0. 0368 

1. 286 

0. 568 

0.  0320 

0. 0400 

1.251 

0.  0229 

0.  0293 

1.  282 

0.553 

0.  0306 

0.0418 

1. 368 

0. 0218 

0. 0303 

1.387 

0.410 

0. 0195 

0. 0247 

1. 264  1 

0.  0139 

0.  01711 

1.  228 

0.  286 

0. 0127 

0. 0155 

1.216 

0.  00910 

0. 01074 

1. 179 

0. 171 

0.  00645 

0. 00767 

1. 190 

0.  00460 

0.00514 

1. 117 

0. 107 

0. 00349 

0. 00448 

1.  285 

0.  00249 

0. 00293 

1. 176 

0.052 

0. 00149 

0.  00204 

1. 366 

0.  00107 

0.  001308 

1. 223 

0.030 

0. 000734 

0. 00109 

1.486 

0. 000524 

0. 00078 

1. 488 

1-in.  Rectangular^ 

3.  D*(I)-! 

a.  17 

a  *  S3.  SS 

0.  808 

0.  0607 

0. 0744 

1.226 

0.670 

0.  0420 

0.0526 

1.252 

4.  0*(1)  =  3#.»» 

o  >3S.Z3 

0. 556 

0. 0307 

0. 0373 

1.  215 

^  =  0.tTC»2 

0.  465 

0. 0229 

0.  0295 

1.  288 

5.  d"(1)  •  55.  K 

0. 340 

0. 0151 

0. 0187 

1.238 

^  =  o.im 

0. 168 

0. 00602 

0. 00765 

1.  271 

0.  071 

0. 00182 

0. 00345 

1.896 

0.0385 

0.000984 

0.00185 

1.  880 

I 
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TABLE  17 


0.  SS^B^E  +  0.  088B^(1  -  S^)  :  EXPERIMENTAL  AND  THEORETICAL  VALUES 


X^(psf) 

or 

exp 

^theor. 

^'exp 

Analytic 

Visual 

.  ®theor. 

18.6 

0. 1255 

0.  201 

0.624 

37.2 

0. 1252 

0. 177 

0.  707 

37. 2  (rect. ) 

0. 1112 

0.119 

0.934 

55.9 

0. 1078 

0. 113 

0. 129 

0.836 

74.5 

0.  0769 

0.  080 

0.  091 

0.845 

Because  of  the  danger  of  direct  beam  contribution,  the  highest  detector  was  an 
inch  below  the  surface,  subtending  a  solid  angle  =  0. 917.  This  is  sufficiently 
close  to  0)  ,  =  1  to  permit  extrapolation  of  the  dose  rate  to  that  value  by  visual  ex- 
tension  of  the  D  (cj)  vs  oj  ciurve  or  by  analytic  curve  fitting  representing  D  by  a 
quadratic  in  u  determined  by  the  method  of  least  squares.  These  curves  are  com¬ 
pared  with  the  data  in  Table  18.  The  results  of  the  two  methods  are  compared  vvith 
the  theoretical  value  in  Table  17.  A  new  set  of  calculated  values  of  D*  is  obtained 
that  at  least  converges  toward  the  experimental  value  at  cu  =  1,  where  agreement 
would  theoretically  be  inevitable,  though  the  discrepancy  at  greater  depths  is  now 
larger.  Values  of  f(d )  for  each  thickness  are  found  to  form  a  fairly  smooth  family 
of  curves  (Figure  9).  A  possibly  significant  factor  that  has  been  neglected  is  the 
variation  of  S  (d,c>))  and  S(d)  with  depth  or  detector  distance  from  the  standard  height 

XL 

of  3  ft.  An  attempt  to  estimate  this  has  been  included  in  Chapter  6. 

General  agreement  has,  therefore,  been  found  between  the  results  of  this  experi¬ 
ment  and  calculations  made  using  the  method  and  curves  of  the  Engineering  Manual . 
However,  significant  discrepancies  exist  that  can  be  accounted  for  theoretically  by 
assuming  that  one  or  all  of  the  factors  B^,  S^,  and  E  are  incorrect  and  that  there 
is  an  additional  factor,  f(d)  that  varies  with  detector  position  and  wall  thickness. 
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0.410  16.71  16.73  16.06  16.07  I  12.34  12.26  8.14  8.50 


CHAPTER  6 


S^(w):  DERIVATION  OF  AN  EXPERIMENTAL  CURVE 

The  infinite  field  dose  rate  has  been  found  for  positions  below  ground  in  a  steel 
cylinder  and  behind  a  rectangular  wall.  These  values  have  been  compared  with 
those  predicted  by  the  method  of  the  Engineering  Manual  using  the  expression 

5^  =  o 

The  factors  B  I  0. 5  S  E  +  0. 088(1  -  S  )  are  represented  by  a  and  can  be  deter- 

IV  I  TV  I 

00 

mined  either  by  using  curves  given  in  the  Manual  or  by  extrapolating  D  to  its  value 
at  o)  =  1  (see  p.  45). 

Since  the  object  of  this  experiment  is  to  evaluate  the  calculated  curve  G  (u)  or 

s 

the  equivalent  S  (u),  it  is  desirable  to  eliminate  the  additional  variables  B  ,  S  , 
and  E  by  following  the  second  course  mentioned  above.  This  experimental  value  of 
or  can  then  be  used  to  find  an  experimental  set  of  values  of  S  (u)  through  the  equation 

a 


^<"u'  '"u> 

Sa<"u'  >  =  -^-5 - ^ 

if  S  (ci)  )  can  somehow  be  evaluated, 
a  u 

For  all  values  of  encountered  in  the  experiment,  is  small  and  so  is  the 
correction  term  S  (u  )  (corresponding  to  the  difference  between  the  experimental 

A  U 

structure  and  an  ideal,  infinite  one).  Further,  this  range  of  S^{u)  can  be  closely 
approximated  by  a  straight  line;  i.  e. ,  S  (u^)  =  Thus  a  simpler  formula  can 
be  found  with  negligible  introduction  of  error 

D"*"^  '“u* 


SO 


I 


t 


The  problem  is  now  reduced  to  evaluating  b.  If  we  return  to  the  original  form 


t 


"u>  =  “ 


=  a 


•"a] 


and  differentiate,  we  find 


dD 


«o 


-  Qfb 


(19) 


(20) 


since  S  (cj  , )  is  independent  of  oj  .  From  the  variation  in  dose  rate,dD‘**/da)  can 
^  u  u  u 

be  found  with  a  change  in  the  height  of  the  cylinder,  corresponding  to  a  change  in 
0)^.  At  a  given  detector  position,  the  height  of  the  cylinder  was  increased  in  steps 
from  1  to  7  ft,  and  measurements  were  made  with  the  source  at  three  radial  posi¬ 
tions,  using  the  1/2-in.  cylinder.  By  dividing  the  change  in  dose  by  the  accompany¬ 
ing  change  in  solid  angle ,  values  for  dD(r)/d(u^  could  be  found  for  each  radial  position. 
An  approximation  of  do'^/do)^  could  be  found  by  adding  the  three  dD(r)/d(i)^  •  Ar  and 
a  fourth  term  representing  the  field  from  45  ft  to  infinity  (Table  19).  This  is  obvi¬ 
ously  a  crude  approach  based  on  marginal  data,  but  it  is  justified  by  the  noncriti¬ 
cality  of  the  exact  value  of  b  and  by  the  consistency  of  the  value  calculated  with  the 

derived  curve  of  S  (w). 

a 

The  theoretical  curve  of  S  (cu)  used  is  taken  from  p.  120  of  Spencer's  monograph 
and  is  based  on  a  constant  detector  height  of  3  ft.  Since  the  experimental  variatimi 
of  0)  was  achieved  by  changing  the  detector  height  in  a  structure  of  c<Hi8tant  dimen¬ 
sions  ,  a  correction  must  be  made  to  relate  the  experimental  data  to  that  which  would 
be  found  at  the  standard  height. 

Using  Spencer's  definition,  we  find 


-1+  w 

^a^“^  “  S^  ^  d(co8  0)  i(d, 
-1 


cos  G). 


B  U 


I 


C 
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The  total  air-scatter  received  from  the  upper  hemisphere  by  an  isotropic  detector 
at  height  d  is  represented  by  S(d)  and  will  decrease  with  height  symmetrically  about 
the  source  plane.  The  integral  term  gives  the  fraction  of  S(d)  scattered  into  the 
solid  angle  of  the  upper  limit  and,  of  course,  equals  S(d)  for  oj  =  1.  Since  detector 
positions  below  ground  receive  higher  dose  rates  than  positions  subtending  equal 
solid  angle  above  the  source  plane,  as  on  p.  45,  it  can  be  assumed  that  Spencer’s 
curve  of  i(d,  cos  0)  will  continue  to  increase  for  values  of  d  less  than  3  ft.  Because 
the  value  of  I  seems  to  increase  with  height  in  the  same  proportion  for  all  values  of 
cos  ©  less  than  0,  the  assumption  was  made  that  the  int^ral  over  t  will  increase  by 
the  same  amoimt.  By  plotting  £(d  -  1)  vs  d  as  given  in  Spencer's  curve  (Figure  10), 
extrapolated  values  of  £  for  d  less  than  3  ft  could  be  found  and  the  corresponding 
values  of 

-1+  03 

y  d(cos  0)  £(d,  cos  0) 

-1 

as  shown  in  Table  20.  Again,  the  inaccuracies  of  the  approach  are  offset  by  the 
smallness  of  the  correction  and  by  the  fact  that  any  correction  will  be  an  improvement. 


Figure  10.  £(d,  cos  0  =  -1)  vs  d 


TABLE  20 

S(d)  and  i(d,  cos  0):  CORRECTION  FACTORS  TO  FIND  D*  ( 3',  cos  0) 

s 


=l(i) 


a 


d 

(in.) 

S(d) 

X(d) 

S(d) 

X(3') 

Correction 

S(3^) 

X(d) 

36 

mgm 

mm 

1.00 

— 

-1 

bb 

1.043 

1.  004 

-2 

0.429 

1. 043 

1.  002 

-2. 75 

0.  088 

0. 4295 

1.043 

0.959 

1.  001 

-3.0 

0. 088 

0. 4295 

1. 043 

0.  959 

1.  001 

-4.25 

0.088 

0.43 

1.043 

0.958 

0.  999 

-4.5 

0.  088 

0.43 

1.  043 

0.958 

0.999 

-5.75 

0.  088 

0. 432 

1.043 

0.  954 

0.  995 

-6.0 

0.  088 

0.432 

1.  043 

0.  954 

0.995 

-8.75 

0.  088 

0.4325 

1.043 

0. 953 

0.993 

-11. 75 

0.088 

0. 434 

1.043 

0.  949 

0.990 

-17.  75 

0. 0863 

0.437 

1.  023 

0.  943 

0. 963 

-23. 75 

0.  0850 

0.44 

1.  007 

0. 936 

0.942 

-35. 75 

0.0844 

0.448 

0. 920 

0. 920 

-47. 75 

0.  0835 

0. 454 

0.989 

0.  907 

0.  898 

COMPARISON  WITH  THEORY 

Once  the  experimental  values  D"’(w,d)  have  been  changed  to  D*’(w,3),  values  of 

S  (cd)  can  be  found  for  each  wall  thickness  (Table  21).  The  values  derived  from 
SL 

these  data  are  plotted  and  are  compared  with  the  calculated  curve  in  Figure  11. 
Curves  through  these  data  points  are  presented  in  Figure  12.  Experiment  and  theory 
come  into  closer  agreement  as  wall  thickness  increases,  indicating  that  the  major 
cause  of  the  discrepancy  lies  in  the  assumption  of  multiple  scattering  in  the  inter¬ 
vening  wall  to  give  essentially  isc^ropic  scattered  radiation.  This  assumption  was 
made^  to  permit  the  use  of  the  solid  angle  subtended  by  the  wall  at  the  detector  as 
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TABLE  21 

S  (w):  EXPERIMENTAL  AND  CALCULATED  VALUES 
a' 


03 

lin. 

1-in.  Rect. 

1-1/2  in. 

2  in. 

Calculated 

0.917 

WBM 

MM 

0.  836 

mm 

mm 

0.  808 

0.680 

■■ 

mgm 

0.776 

0. 643 

0.677 

0.674 

■PH 

0.  757 

0. 620 

0.615 

mm 

0.639 

mgm 

0.670 

0. 482 

0. 

0.502 

0. 543 

0.492 

0. 504 

mgm 

0.649 

0.502 

0. 483 

0.496 

0. 488 

0. 568 

0. 406 

0. 418 

0.379 

0. 390 

mgm 

0. 556 

0. 343 

0. 553 

0. 398 

0. 394 

0.395 

0. 400 

■■■ 

0. 465 

0.  273 

0.225 

0.410 

0. 275 

0. 265 

0.236 

0.  230 

0. 185 

0.  340 

0. 174 

0.145 

0.286 

0. 194 

0. 178 

0.150 

0. 146 

0. 122 

0. 171 

0. 105 

0. 0945 

0. 0753 

0. 0712 

0.0628 

0. 168 

0. 0745 

0.0620 

0. 107 

0.  0648 

0. 0549 

0.  0418 

0.  0350 

0. 0710 

0.  0343 

0. 0218 

0. 052 

0.  0317 

0.  0262 

0. 0219 

0. 0201 

0. 0159 

0.  0385 

0. 0195 

0. 0112 

0.030 

0.0183 

0.0151 

0.  0127 

0.  0127 

0. 0084 

the  sole  criterion  of  the  radiation  scattered  to  the  detector,  whether  from  above  or 
below  the  detector  plane.  Obviously,  this  assumpticm  will  be  more  reascmable  in  a 
structure  with  thicker  walls  in  which  multiple  scattering  will  more  effectively  con¬ 
ceal  the  different  angular  distributions  of  the  radiation  falling  on  different  sections 
of  the  wall. 


Figure  9.  in  which  f(d)  is  plotted  for  each  value  of  X^,  gives  an  indicati<m  of 

the  variation  of  the  scattering  properties  of  the  waU  with  detector  positions.  This 

is  essentially  the  same  as  a  ratio  of  the  experimmital  and  calculated  values  of  S  (u) 

ft 

given  in  Table  21.  A  series  of  curves  for  S  (o))  has  thus  been  found  that,  in  tiie 

ft 

limit  of  great  thickness,  approaches  that  calculated  Iqr  Spencer.  They  differ  sub- 
stantially  for  the  smallest  solid  an^es  and  thinnest  walls,  though  this  discrepancy 
is  somewhat  o^et  in  dose  calculations  by  a  discrepancy  in  a  in  the  opposite  direc¬ 
tion.  Since  experiment  predicts  a  higher  dose  rate  than  theory,  substantiating  ex¬ 
periments  would  be  desirable. 

This  experiment  indicates  that  the  calcnlational  methods  of  the  Engineering 
Manual  are  entirely  adequate  for  the  evaluation  of  most  shelter  situations.  However, 
there  does  seem  to  be  room  for  a  considerable  refinement  of  theory  to  interpret  the 
discrepancies  encountered. 
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